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We apply a particle-particle RPA model to study the properties of the two-neutron valence 
wave function in nuclei 14 C, 12 Be, n Li and 14 Be. The RPA model takes account of two- 
body correlations in the cores so that it gives a better description of energies and amplitudes 
than models which assume a neutron closed shell (or subshell) core. With a Gogny neutron- 
neutron effective interaction or with the equivalent density dependent delta force we are able to 
reproduce the two-neutron separation energies in these nuclei and in the corresponding cores, 
except for 9 Li. These calculations suggest the same 2s-lp;i/2 shells inversion in 12 Be- 13 Be than 
in 11 Be. 


I. INTRODUCTION 

In an earlier work Q, we calculated the properties of the two-neutron valence pair in the nuclei n Li, 12 Be, 
14 C using a Tamm-Dancoff model and assuming the core to have a closed IP 3/2 neutron shell. Important 
ingredients in the calculations are the assumed pairing interaction and the role of the P 1 / 2 -S 1/2 shells inversion 
that is visible in the spectrum of 11 Be. 

While the model gave a coherent and reasonable description of these nuclei it fails to describe 6 He (taking 
now an alpha-particle core) giving a two-neutron separation energy of several MeV instead of the experimental 
value of 0.97 MeV j|. This has already been found by Esbensen et al. j|. When one compares 6 He described 
as an alpha-particle + two neutrons to 14 C for example, one sees immediately an important difference between 
the two systems: for 6 He the core is well described in an Hartree-Fock model as a pure (ls) 2 (ls ) 2 configuration 
||] while we have long known Q that in 14 C the core of 12 C is a mixture of states with neutrons in the IP 3/2 
or lpi /2 states, what is not taken into account in Tamm-Dancoff models. 

In this work, we will include the core correlations using the two-particle RPA theory and also examine more 
broadly the dependence on the assumed effective two-neutron force. Zero range residual forces are very simple 
to use but are quite arbitrary and binding energies are not sufficient to fix uniquely the force. In a recent work 
Garrido et al.jG) have searched for a zero range density dependent force equivalent to the finite range Gogny 
interaction (see also [||). The authors fit the force in order to reproduce the gap calculated in nuclear 
matter with the Gogny force which has good pairing properties in finite nuclei. A fit to the whole domain of 
kf values determines unambigously the parameters of the force and tells what is the cut-off on neutron energy 
to be used. This last information is very important since a zero range interaction has no natural cut-off. The 
density independent part of the force reproduces the low energy properties of a free neutron-neutron system 
so that their force is equivalent to a realistic effective interaction in two-neutron and infinite systems. We will 
use this force in our problem and will discuss the results compared to the zero range pairing forces used in 
Tamm-Dancoff models. 

In section II we briefly report on the results obtained in a pairing or Tamm-Dancoff model with three effective 
neutron-neutron interactions. In section III we recall the properties and equations of the particle-particle RPA 
model. The results of this model are presented and discussed in section IV for 14 C- 10 C, 12 Be- 8 Be and 11 Li- 7 Li 
and in section V for 14 Be- 10 Be with a discussion on the 13 Be and n Be spectra. At the end, section VI is devoted 
to our conclusions. 


II. HAMILTONIAN AND EFFECTIVE INTERACTIONS 


We first make a pairing, or Tamm-Dancoff, approximation to describe core + two neutron systems. Assuming 
an inert and closed sub-shell core for neutrons we diagonalise the two-neutron hamiltonian in a two-neutron 
subspace built on non occupied neutron states in the core: 
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h nc (i) — + Kic(*) 

2 /r 


2/, xlrf/W 


KcW = —V NZ ( /(f) - 0.44r-5(l.s)- ^ 


r - R 0 


16 a 2 a r , 


df(r) 

dr 


f(r) = ^1 + exp( 
Vnz = Uq — U-, 


a 

N - Z 


(3) 

(4) 

(5) 

( 6 ) 


fi is the reduced mass equal to 
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the core; Rq = tqA c with 7'o=1.27 fm, a=0.75 fm; the strengths Uq and U T are taken the same as in our 
previous papers |I|.|To|. The last term of potential, eq.(4), simulates particle-phonon couplings contribution to 
the one-body potential Jn[]. The strengths a n are fitted in each nucleus to reproduce the experimental lpi/ 2 , 
2s, ld 5/2 single neutron energies for 14 C and 12 Be and in order to get the measured two-neutron separation 
energy in 11 Li and 14 Be for lpi /2 and 2s states which are not experimentally well known . Their numerical 
values can be found in refs.[l] and [10]. For higher neutron states the particle-phonon couplings are weak and 
we take a n =0. We use a discretisation of the continuum states with a radial box of 20 fm and orthonormalise 
the wave functions by the Schmidt method. In our previous papers the two-body term P 1 .P 2 was neglected but 
as the approximation was the same for all nuclei: 14 C, 12 Be, n Li and 14 Be and the effective neutron-neutron 
interaction fitted to describe the properties of 14 C and 12 Be, then used in the other nuclei, the effect of such 
an approximation was minimised. We have checked that adding this term changes slightly the strength of the 
effective interaction (by few %) but gives the same agreements and predictions than previously. 

We used in refs.[1,10] a zero range density dependent neutron-neutron interaction given by: 


V„ n (l,2) = -Vo 1 -x[p c ( 


r ri +r 2 , 


)/Po] p <5(Tl - r 2 ) 
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where p c is the core density and po=0-16fm -3 . The above fit gives now Vo=880 MeV.fm 3 , x=0.93 and p=1.2 
close to what is used in the literature |],|l2j]. With the same force we have calculated 6 He as an alpha-particle 
plus two neutrons. Taking V nc as the neutron-alpha particle potential fitted by Satchler et al. ]l3| to low 
energy 1=1 phase shifts we get a two-neutron separation energy of about 4 MeV in 6 He while it is 0.97 MeV 
experimentally Similar strong binding was found by Esbensen et al. || who modify the n-n interaction 
in order to get the experimental value. However this is not satisfactory for us because we always require that 
our effective interaction leads to an overall agreement for all systems since the density dependent term should 
take account of the dependence of the effective interaction upon the nucleus. We have looked for different 
possible sources of such unexpected disagreement: choice of neutron-alpha particle interaction, discretisation 
of the continuum which will be discussed in a forthcoming paper and effect of the cut-off on neutron energy. 
However none of these effects is responsible for such a bad result. 

The force of eq.(7) has three parameters which are not independently determined from a fit of energy spectra. 
In a recent paper Garrido et al. have fitted the three parameters of V nn in order to get the same gap, A (kp), 
in nuclear matter than obtained with a Gogny finite range effective interaction |^,||. They have shown that to 
get agreement over all the domain of kp they have to take p=0.47, x=0.45 with a cut-off, e®, on the neutron 
energy of 50-60 MeV. The cut-off energy determines the strength Vo if one assumes that the density independent 
part of the interaction should reproduce the properties of a free neutron-neutron system. It gives the following 
relation between Vq, e° and the neutron-neutron scattering length a nn 0 : 
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a nn is experimentally -18.5 fm then very large. If we replace it by -cso and take e?=60 MeV this relation gives 
V 0 = 480 MeV.fm 3 as used in ref. [6], 

For free neutrons, the neutron energy is only kinetic energy while in finite nuclei the neutrons are in the 
potential V nc . Therefore when using eqs.(8-9) the cut-off energy has to be calculated from the bottom of the 
single-particle well and e 3 replaced by : 


2 











e c + Vnz 


( 10 ) 


e 


o 

c 


where Vnz is the depth of our one-body potential, eq.(6), and e c the cut-off energy for a neutron described 
by the hamiltonian of eq.(3). In our calculations we take neutron states up to an energy e c ~ 10 MeV for all 
nuclei. Vnz depends on proton and neutron numbers in the core and varies between about 40 MeV for n Li 
and 50 MeV for 14 C , giving an equivalent cut-off in free two-neutron system or nuclear matter of 50-60 MeV as 
required by the fit of ref. [6]. Note that taking e c the same for all nuclei implies that Vq will be slightly different 
for 11 Li, 12 Be and 14 C. 

This new zero range effective interaction has very different parameters compared to the usual ones. We have 
made the same calculations as previously with this force for our systems 14 C, 12 Be and 6 He. n Li results are 
not reported in the table but with neutron energies close to experimental values it is not bound. The results 
are reported in Table I where they are compared to the measured two-neutron separation energy, [^|. We 
see two interesting facts: the calculated S^n in 6 He is now 1.02 MeV close to the experimental value but it is 
much too low in the other nuclei. Since this zero range force is constructed to reproduce the Gogny force in 
nuclear matter and in a free two-neutron system we thought that this equivalence could fail in our finite nuclei 
which are quite far from both systems. 

Then we have made the same Tamm-Dancoff calculation with the genuine Gogny force, D1 or D1S ]t]||. The 
Gogny forces are the sum of central density dependent and spin-orbit terms. For two neutrons coupled to 0 + the 
spin-orbit term gives very small contribution to pairing matrix elements and can be neglected 0 - Therefore 
we have made our calculations, keeping the central term only. The results are very close for the two forces and 
we present them for the force D1S only. The results with the Gogny force are presented in Table I. 


Force 


i2 Be 

b He 

Gogny 

11.8 

2.23 

0.94 

<5 force 

11.7 

1.95 

1.02 

exp. 

13.12 

3.67 

0.97 


TABLE I. Two-neutron separation energies in MeV calculated in Tamm-Dancoff model with the Gogny DS1 
force and a zaro range force of eq.(7) with Vo=480 MeV.fm 3 , p=0.47 and x=0.45 


We see that they are the same as obtained with its zero range substitute what tells us that: 

1 - the equivalence between the two forces, with zero or finite range, shown for a free system and nuclear 
matter holds in all considered finite nuclei. We give the results for Lo=480 MeV.fm 3 , the same for all nuclei, 
but as already mentionned Vq should be larger in 12 Be than in 14 C following relations (6), (8) and (10) what 
will improve the agreement between the two series of results. 

2 - the range of the effective force is not responsible for the unability of the nuclear model to describe simul¬ 
taneously 6 He and the p-shell nuclei. 

The agreement for 6 He on one side, the disagreement for the other nuclei on the other side seem at first 
sight a contradiction of the model. However it can be understood when we go back to the first assumption of 
the Tamm-Dancoff model that the core is inert with neutrons filling the lowest shells. For 6 He the core, an 
alpha-particle, is strongly bound and very well described in an Hartree-Fock model with neutrons and protons 
filling the Is shell Q]. It means that the assumption of a closed shell core (implicitly assumed when we put 
our extra- neutrons on the unoccupied shells with probability one) is valid in this case. However for the other 
nuclei we know that the cores are not good closed shell nuclei. It is known for long in the case of 12 C from the 
work of Cohen-Kurath For 12 Be shell model calculations also show a large deviation from a pure state as 
it is required by experiments [fi"5| Il8|. It is also obvious in our calculations for 12 Be (ll or in three-body Fadeev 


model calculations of Thompson et al. |lTJ where 12 Be is described as a core of lb Be plus two neutrons and 
found as a large mixture of (2s) 2 , (lpi/ 2 ) 2 and (ld 5 / 2 ) 2 two-neutron components while it is considered as a pure 
(2s) 2 or (lpi/ 2) 2 state when one studies 14 Be. 

A model which takes account of two-body correlations in the ground state of the core is the particle-particle 
RPA model more generally called pair vibration model [20 2^]. We briefly recall what this model is for two 
neutrons outside a core which, in an independent particle model, would have closed shells (or subshells) for 
neutrons. 


III. PARTICLE-PARTICLE RPA MODEL 

Particle-particle RPA relies on the expansion of the two-particle (two-hole) Green’s function in terms of ladder 
diagrams with upward and backward going diagrams as shown in Fig.l. Note that a summation over upward 
going diagrams only leads to the Tamm-Dancoff approximation of the previous section. From the related 
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approximate integral equation satisfied by this RPA two-particle Green’s function, a system of equations is 
derived which describes simultaneously the A c +2 and A c -2 systems where A c characterises the core which in 
an independent particle model would be described as a closed shell (or sub-shell) nucleus. 



FIG.l: RPA diagrammatic expansion of the two-particle Green’s function. 


We here apply this method to our problem of two neutrons outside a core with N c neutrons assuming that 
the protons are not disturbed when one adds or subtracts two neutrons. 

Let’s define a=(ai,a2), b=(bi,b2)... two-neutron configurations with the neutrons in states ai,a2... unoccupied 
in the Hartree Fock core ground state, a,/3... two-neutron configurations with the neutrons in occupied states 


ai,o;2--- and two-neutron amplitudes or spectroscopic factors as: 

X a (J, M) = < N C + 2\A+{J,M)\N C ,0 > (11) 

X a (J,M) = < N c + 2|A+ (J, M)\N C , 0 > (12) 

where the pair operators are given by: 

A t( J , M ) = Y (iai)ia2> TO o 1 ,m 02 |J,M)a+a+ with oq < a 2 (13) 

m ai ,rria 2 

A i(j,M)= y (ja 1 Ja 2 , m a 1 ,rn a 2 \J,M)a+ 1 (ia2 with ai < a 2 (14) 

m ai ,m a2 


ail i s th e creation operator of a neutron in state ji, rrn, |-/V C + 2 > and \N C , 0 > are respectively the RPA wave 
functions of the N c + 2 nucleus in excited or ground state and of the core in its ground state. We see immediately 
that the amplitudes X a are non zero only if the core ground state has 2p-2h,4p-4h,... components. Assuming 
that all X a are zero, one gets back to the Tannn-Dancoff approximation. 

In the same way we define two-neutron hole amplitudes for the N c — 2 nuclei from A a and A a , the annihilation 
operators for a pair which are hermitian conjugates of A+ and A+ defined by eqs.(11-12). They are: 

Y a (J, M) =< N c - 2\A a (JM)\N c , 0 > (15) 

Y a {J,M) =< N c -2\A a (J,M)\N c ,0> (16) 

For the N c -2 nuclei Y a are non zero only if the core ground state is not a pure HF state but has 2p-2h, 4p- 
4h... components. Then these anomalous components X a and Y a give a measure of the two-body correlations 
in the cores. For a given spin and parity the RPA amplitudes X and Y satisfy the same system of equations: 

(E - e a )x a Y K a \ V nn\b > X b - Y < a \ V nn\P > Xff = 0 (17) 

b 0 

(E - e a )x a + Y < a \Vnn\b > Xb + Y, < a \^nn\f3 > X0 = 0 (18) 

b 0 

where x are the amplitudes X or Y as explained below. The two-body matrix elements are antisymmetrised and 
e a , e Q the sum of unperturbed energies of the two neutrons in configurations a, a respectively. The one-neutron 
states are described by the one-body effective potential of eq.(4) what means that we expand our two-body 
Green’s function of Fig.l in terms of one-body propagators which include particle-phonon couplings as given in 

Fig-2. 
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FIGURE 2. One-particle Green’s function with inclusion of particle-phonon coupling diagrams. 


If we take N configurations a,b.. and M configurations a , the equations (17,18) have N+M solutions. N 
solutions correspond to the N c + 2 nucleus with eigenvalues E n j and amplitudes where n labels the state 
such that: 


E nJ = E nJ (N c + 2)-E 0 (N c ) (19) 

X<?HJ)=xW , A(20) 

the amplitudes satisfy the following orthonormalisation relation: 

£ x™ Xi n>) £ £ n) xW = <W (21) 

a a. 

The M other solutions (labelled by the index nr) are eigenstates of the N c — 2 nucleus with eigenvalues E m j 
and amplitudes Y given by: 


£ m j = ( E m j(N c — 2) — E 0 (N C )) = —E mJ (22) 

£"£)=*£ ,Y^\j)=x^ (23) 

with the orthonormalisation condition: 

£ yw y ,£> - £ yj ra '> = -<w (24) 

a a 

We see from eqs.(19) and (22) that the lowest energies E 00 and £ 0 o give the two-neutron separation energy in 
the N c + 2 and N c nuclei respectively with the relations: 


S 2n (N c + 2) = -E 00 ( 25) 

S 2n (N c ) = £oo (26) 

Then, contrary to particle-hole RPA, all solutions of eqs.(17,18) correspond to physical states. The separation 
between N and M states follows from the relative importance of x a and x a amplitudes: for the N c + 2 nucleus 
the amplitudes x a are larger than x a and inversely for N c — 2 nucleus. 

This model is now applied to the p-shell nuclei. As noticed already 6 He described as an a particle plus two 
neutrons is very likely well reproduced by a Tamm-Dancoff model since the core is a good Hartree-Fock system. 
On the other hand, the Ay-2 nucleus (di-proton) is not a bound system and the RPA model is not reliable in 
this case so that we discuss the p-shell nuclei only. 


IV. RESULTS FOR N=8 NUCLEI: 14 C, 12 Be AND n Li 

We make the calculation for the two effective interactions, the Gogny finite range interaction and the zero 
range interaction of eq.(7) with the parameters p and x fitted by Garrido et al. and a strength V$ fitted to give 
the same RPA energies as the Gogny force. Afterwards we shall compare with the strength given by eqs.(8-10). 
The two-neutron subspace (a,b...) is the same as in section II with a normal ordering of lpi/ 2 - 2 s shells in 
C-isotopes but an inversion of these two shells in Be and Li. For the occupied neutron states (a, /?...) we take 
the two shells lsi /2 and IP 3/2 ■ The IP 3/2 neutron energy is taken as the known one-neutron separation energy 
for the corresponding core ( 12 C or 10 Be or 9 Li ) Q . For the Is shell we have no experimental information and 
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we take it as a parameter by changing the depth of our Saxon-Woods potential. This very deep state has very 
little effect on the N c +2 states but a non negligible effect on the energy of the N c -2 ground states. Then we 
take the Is energy in order to get close agreement, if possible, with the experimental two-neutron separation 
energy in the core nucleus. This procedure gives e(ls)=-32, -32 and -28 for 12 C, 10 Be and 9 Li respectively. 



14 C 1«C 

12 Be iu Be 

xl Li y Li 

Gogny 
6 force 
exp. 

12.9 32. 
12.9 32.5 
13.12 31.8 

3.69 8.2 
3.6 8.52 
3.67 8.48 

0.37 3.8 

0.34 3.76 

0.34 ±0.05 6.1 


TABLE II. Same as Table I for RPA model and Vo=500, 510 and 560 MeV.fm 3 for 14 C- 12 C, 12 Be- 10 Be and 

11 Li- 9 Li respectively 


In Table II are presented S 2 „, the two-neutron separation energies in N c +2 ( 14 C, 12 Be and 11 Li) and N c 
( 12 C, 10 Be and 9 Li) nuclei given by the lowest eigenvalues Ego and £oo respectively, following relations (25) and 
(26). The results are given for the Gogny force and the equivalent zero range force. In Table III we give the 
amplitudes X a and X a for N c + 2 nuclei and Y a and Y a for N c — 2 nuclei corresponding to the ground states 
and obtained with the Gogny interaction. The zero range interaction gives nearly identical eigenvectors. 

For 14 C- 12 C and 12 Be- 10 Be, S‘ 2 n is close to the experimental values while we have no free parameters, apart 
from the ls-energy which turns out from the fitting procedure to be very close to a typical Hartree-Fock energy. 
We see that the amplitudes X a for 14 C and 12 Be are for a = (I.P3/2) 2 rather large and larger in 12 Be than in 
14 C. This large value means that in the ground state of the cores, 12 C and 10 Be, there are large components 
of 2nh-2np states with at least two holes in the lp 3 / 2 shell. This is in complete agreement with the shell 
model results for 12 C || and for 12 Be O [li|. Furthermore we see that, if in 14 C the amplitude for the two 
neutrons in a lpi/2 state is very large, i~ 2 Be we have comparable amplitudes for (2s) 2 , (lpi/2) 2 and (ld.5/2) 2 
configurations. This means that, together with the fact that X a = 0.46, 12 Be cannot reasonably be considered 
as a closed shell nucleus as done in many papers on 14 Be. 

Our RPA equations give simultaneously the amplitudes Y a and Y a of eqs.(15-16) for 10 C and 8 Be. Here the 
Y a are the anomalous amplitudes reflecting again correlations in the cores. We see consistency with the results 
for 14 C and 12 Be. The amplitudes Y„ are larger in 8 Be than in 10 C with a distribution over several two-neutron 
states revealing a very complicated structure of 8 Be. 

The u Li- 7 Li systems present more ambiguity than the previous ones. Now we know from break-up ex¬ 
periments that the lowest neutron resonance in 10 Li is an s-state at 0.1-0.2 MeV |23|as required by several 
calculations sum and that the next state is a p!/ 2 resonance at 0.54± 0.06 MeV ||26|. 



X a (Y a ) 

X a (Y a ) 


(2s) 2 

O 1 / 2) 2 

(lpi/2,2pi/ 2 ) (ld 5 / 2 ) 2 

(IP 3 / 2) 2 

(Is) 2 

~ ra C r 

-0.12 

0.96 

-0.05 

-0.28 

0.19 

-0.09 

10 C 

-0.05 

0.15 

-0.02 

-0.14 

0.98 

-0.32 

TT Be 

-0.48 

0.76 

-0.29 

-0.44 

0.60 

-0.09 

s Be 

-0.16 

0.34 

-0.17 

-0.30 

1.18 

-0.12 

n Li 

0.66 

-0.56 

0.48 

0.04 

-0.45 

0.07 

7 Li 

-0.12 

0.22 

-0.27 

0.02 

1.12 

-0.11 


TABLE III. The most important RPA amplitudes (X a , X a for N c + 2 nuclei, Y a , Y a for N c — 2 nuclei) 


The results, presented in Tables II and III, have been obtained with e(ls) = —28 MeV, e(2s) = 0.19MeV 
and e(p!/ 2 ) = 0.6 MeV close to measurements for the two last ones and reasonable for the first one. We first 
see that Si n in 11 Li is in good agreement with the experimental value, 0.34±0.05 MeV j2?j]. The anomalous 
amplitude X a for ai=(lp 3 / 2 ) 2 is as large as it was in 12 Be revealing strong deviation from closed shell in 9 Li. 
We also see that the amplitudes for (lpi/ 2 ) 2 and (2s) 2 configurations are similar, as required by the analysis of 
break-up experiments [^8| . However the Si n value for 9 Li is very far from the measured value, 3.85 instead of 
6.1 MeV. This disagreement, found only in Li, cannot be improved by a reasonable change of e(ls) which is the 
only parameter of the calculation. It could come from the single last proton which is on the same shell as the 
active neutrons in 'Li and 9 Li but we did not find any simple way to evaluate such an effect. 

In the case of Li isotopes, there is an ambiguity coming from higher neutron states. Experimentally one has 
not seen the ds/ 2 resonance, therefore we have calculated it with the bare Woods-Saxon potential even though 
we know that in n Be one has to add a surface term which lowers the d 5 / 2 energy. However the effect of such 
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a term should not affect the results since, as seen in Table III, the amplitude for the (lds/ 2 ) 2 configuration 
is very small. Also, because 10 Li is unbound, all neutron states are in the continuum with the result that 
continuum non-resonant states play a more important role in the determination of the 11 Li energy. This effect 
is particularly important for the discretised 2p 1/2 state which is taken as an eigenstate of the bare Woods-Saxon 
potential. Because the modification of the potential by particle-vibration couplings is very large for the lpi /2 
resonance, this 2 p !/2 state comes at an energy close to the lpi/ 2 , giving a large overlap between the two states. 
We know that particle-vibration couplings are more important for states close to the Fermi surface and can be 
neglected for higher states in usual nuclei. However in 10 Li these couplings are very strong and very likely still 
efficient for the 2p]y2 state. Then as a test we have made the calculation with the same surface contribution to 
the one body potential for lpi /2 and 2p!/ 2 states, taking the bare potential for the other p-states which anyway 
have much higher energies and therefore have a weaker effect on the ground state energy and wave function. 
The amplitude X a for a=(lpi/2,2p!/ 2 ) is weaker, S^n slightly smaller. To recover the same value one has to 
take e(lp!/2)=0.53 MeV, then at the experimental energy, but these differences are not significant. 

Going back to Table II we see that the strengths Vq of the zero range interaction fitted to recover the same 
energy than the Gogny interaction are 500, 510 and 560 MeV.fm 3 for 14 C, 12 Be and 11 Li respectively. They 
follow closely the dependence on Vnz given by eqs.(8-10) with a value for 14 C very close to 480 MeV.fm 3 used 
in ref. [6]. We again find the same equivalence between the two forces,with zero or finite range, as derived for 
free neutrons and nuclear matter. It seems from our study that the need of a much stronger zero range force to 
get experimental binding energies in simple pairing models, as reminded in Section II, is in fact a way to take 
implicitly account of two-body correlations in the core which are neglected in the model. 


V. RESULTS FOR 14 Be (N=10). 

The problem of 14 Be is not yet completely clarified, both experimentally and theoretically. Experimentally 
one knows that 13 Be has a d 5/2 resonance at 2.01 MeV above the n+ 12 Be threshold J29|, a lower resonance 
at about 0.8 MeV seen in the reaction 14 C( 11 B, 12 N) 13 Be |io) which has no spin or parity assignment and, 
from a recent experiment | pl| | using fragmentation of 18 0 and detecting neutrons in coincidence with 12 Be, a 
l/2 + resonance below 0.2 MeV which should be the ground state of this unbound nucleus. Theoretical models 
describing 14 Be as two neutrons outside a core of 12 Be, where the neutrons fill the Is, IP 3/2 and lpi /2 shells, 
have either to lower the d 5 / 2 resonance or to assume a bound 2s state therefore to bind 13 Be to get a correct 
binding energy in 14 Be & In a first paper ||lG[] based on a two-neutron Tamm-Dancoff model and with a zero 
range effective interaction fitted on 14 C and ^^Be we found that an inversion of the lpi/ 2 - 2 s shells leads to the 
correct binding energy in 14 Be without modifying the known ds /2 energy nor assuming a bound 13 Be. This 
calculation, as others, assumes a closed shell nucleus of 12 Bc where the neutrons fill the ls-lp 3 / 2 -2s shells, while 
we have seen above that the ground state of 12 Be has amplitudes X a over several two-neutron configurations 
and, because of large X a , has very likely components on more complicated configurations. The particle-particle 
RPA applied to 14 Be takes this into account and will provide a description of the core consistent with the RPA 
amplitudes of Table III. 

We have made the RPA calculation assuming a normal ordering of shells: Is, IP 3 / 2 , lpi/ 2 , 2s,...with a 
Hartree Fock state where the p-shell is filled for neutrons. With e(lp 1 / 2 )=-3.15 MeV, given by the neutron 
separation energy in 12 Be [|2|, a ld 5 /2 state at 2 MeV, the experimental value, and a 2s state at 7 keV, very 
close to the threshold, we get S 2 n = 0.7 MeV much too low compared to the experimental value of 1.34±0.11 
MeV [j2j. For 12 Be we get S2 n =3.26 MeV, also too small. Moreover the amplitudes X a are small indicating 
weak correlations in 12 Be in disagreement with the results of the previous section. Therefore we find again that 
we are not able to describe satisfactorly 14 Be when we assume a normal ordering of shells in 12 Be. 

We now assume an inversion of the two shells 2s and lpi /2 as in n Be. The configurations a are built on 
ls,lp 3 / 2 and 2s states while the configurations a include neutrons in lpi /2 state. The results are summarised 
in Table IV where we give S 2 ra for 14 Be and the RPA amplitudes for 14 Be and 10 Be. They are obtained for a 
Gogny effective interaction with the d 5 /2 resonance at the experimental energy of 2 MeV, a lpi /2 resonance 
at 0.68 MeV, an occupied 2s shell with an energy of -3.15 MeV given by the experimental neutron separation 
energy in 12 Be and a lp 3 /2 state at -5.6 MeV. For the a configurations we have taken the (IP 3 / 2 ) 2 , (Is,2s) and 
(2s) 2 states. We have checked that adding the (Is) 2 states does not change our results. 

We get 3.65 MeV for the two-neutron separation energy in 12 Be. By comparing with the result of Table II 
we see that the values of S 2 ra in 12 Bc found in both calculations are very close showing the coherence of the 
model. Indeed we get 3.65 MeV when 12 Be is considered as the core for the 14 Be- 10 Be systems while it is 3.69 
MeV when it is described as two neutrons outside a core of 10 Be. 

The p 3 /2 energy cannot be deduced directly from experiments. However n Be has a 3/2“ state at 3.9 Mev 
excitation energy |3^], then at 7 MeV when referred to the 12 Be core, with a very small width and has certainly 
a large component of one neutron-hole in 12 Be mixed with 2h-lp components which have higher energies. 
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Therefore a hole energy of 5.6 MeV is not unrealistic. 

According to our discussion for n Li we have done the calculation using the same surface potential for lpi /2 
and 2p!/ 2 states. We have then to lower the lpi /2 energy to 0.56 MeV if we leave all other states to be the 
same , what does not modify qualitatively our conclusions. 



S 2 n(MeV) 

(lPi/ 2) 2 (lpi/ 2 , 2 pi/ 2 ) (ld 5 / 2 ) 2 

(lp 3/2 ) 2 ( 2 s ) 2 

"^Be 

10 Be 

1.30 

-0.46 0.62 0.74 

-0.20 0.31 0.46 

-0.64 0.46 

-0.73 1. 


TABLE IV. RPA energy in MeV for 14 Be and main amplitudes for 14 Be and 10 Be obtained with the D1S 

Gogny force. 


The amplitudes given in Table IV show again strong correlations in 12 Be. Indeed the anomalous amplitudes in 
14 Be are 0.64 and 0.46 for (P 3 / 2) 2 and (2s) 2 respectively what means that in 12 Be components on configurations 
with two holes on these shells are important. Note that the components with two holes in the 2s-shell are 
qualitatively consistent with the amplitudes X a of Table III for a=(lp!/ 2 ) 2 and (ld 5 / 2 ) 2 - However RPA gives 
only overlap of wave functions and a strict and direct comparison between amplitudes derived in the two 
calculations is not possible. To make a direct comparison one has to calculate wave functions what can be made 
using a quasi-boson approximation as was done in the past for particle-hole RPA correlations. 

Our fitted value e(lpi/ 2 )= 0.68 MeV is consistent with an unbound 1/2 - state in 13 Be which would be at 
0.68 MeV above the 12 Be +n threshold, close to the experimental resonance at 0.8 MeV. However the recent 
experiment using fragmentation of 18 0 |3l]] shows a low energy (<0.2 MeV) s-wave strength in 13 Be what, in 
an independent particle model, would mean that the lowest unoccupied shell in 12 Be is an s shell and would 
reject the possibility of inversion. In our model however 13 Be is described as a neutron added to the correlated 
core of 12 Be which is a mixture of many different states, in particular it has large components on configurations 
with a closed IP 3/2 shell plus two neutrons on the 2s or on the lpi /2 shell. Therefore these two components 
will give in 13 Be two different states, a 1/2“ state built on the first one with the last neutron on the P 1/2 shell 
(the 2s-shell is filled) and a l/2 + state built on the second one with the last neutron on the empty 2s shell. In 
a weak coupling model, because of the known inversion in n+ 10 Bc system, one can show that the l/2 + state is 
lower than the 1/2“ state by about 0.32 MeV and at about 0.3 MeV above the 12 Be+n threshold what is in 
qualitative agreement with the experimental state seen recently. Consequently a l/2 + ground state in 13 Be is 
not in contradiction with an inversion of the 2 s-lpi /2 shells. 

There are several other arguments in favor of this inversion in 12 Be- 13 Be. The first one relies on the recent 
measurement of the B(E2) for the 2 + state at 2.1 MeV in 12 Be which is found to be the same as in 10 Be for 
the 2 + state at 3.4 MeV |k|. Because the inversion in n Be is related to the large value of the B(E2) in 10 Be 
[11 |,b 3 there is no reason why the effect should be smaller in 13 Be. Moreover the phonon has a smaller energy 
in^^Be than in 10 Be what is expected to give an enhancement of the coupling Jll|]. Note that this large B(E2) 
was predicted in ref.[10]. Further arguments are found by looking at the n Be spectrum. Indeed if 13 Be can 
be described as a neutron added to a core of 12 Be, n Be can be described as a hole in the same core of 12 Be. 
Therefore the l/2 + ground state of n Be is expected to correspond mainly to a hole in the last occupied shell 
in the Hartree-Fock ground state of 12 Be what implies that this shell should be an s shell, not a Pi/ 2 - The 
first excited states can be obtained as two holes coupled to 0 + plus a neutron on the p x / 2 or the d 5 / 2 shell 
therefore will be a 1/2“ and a higher 5/2 + states respectively, in agreement with the experimental spectrum. 
If we assume that the residual interaction between the two holes and the particle are similar if the particle is 
on a P 1/2 or d 5 / 2 shell we find a difference between the excitation energies of the two states of 1.32 MeV while 
experimentally it is 1.45 MeV. One sees that the inversion is able to give a coherent description of n Be, 13 Be 
and 14 Be . 

Because any model relies on approximations one should be aware of the difficulty to draw a precise scheme. 
For this reason it is desirable to compare results of different models. Our results without inversion of shells 
agree with those of Thompson and Zukhov Jlff l . We now make a comparison with the work of Descouvemont 
|35],[j(|. He has calculated in the generator coordinate model (GCM) 13 Be and 14 Be as 12 Be +n and 12 Be -kn+n 
systems respectively. The core of 12 Be is described by a filled lp-shell for the neutrons while the two protons 
of the lp-shell can couple to different states, ground and excited 0 + , 1 + and 2 + states. Qualitatively it is 
equivalent to our Tamm-Dancoff approach where contribution of core excited states are put in our one-body 
neutron potential and where the neutrons of the core are assumed to fill the lp-shell. The GCM calculations 
lead to a slightly bound l/2 + ground state for 13 Be and a 14 Bc bound by 1.1 MeV. In our Tamm-Dancoff 
approach with the zero range force fitted on 14 C and 12 Be, a 2s neutron state at -90 keV and a filled neutron 
lp-shell (therefore without inversion) we get for 14 Be a binding energy of -0.93 MeV what is close to the GCM 
result. The three models, three-body Faddeev, GCM or simple pairing, lead to similar results but are not able 







to give rise to good agreement for 13 Be and 14 Be when in 12 Be the neutrons are assumed to fill the IP 3 / 2 -IP 1/2 
shells. 

A different work |5?j using a density dependent relativistic mean field model calculates 12 Be- 14 Be assuming 
closed lp and lp-2s neutron shells respectively. It gives too large one- neutron and two-neutron separation 
energies in both systems. This result, together with the assumption of closed shells in both 12 Be and 14 Be, is 
not in agreement with other calculations. 


VI. CONCLUSIONS 


We have first shown that a two-neutron Tamm-Dancoff model with a zero range density dependent neutron- 
neutron interaction fitted on 14 C and 12 Be gives simultaneously good results for n Li and 14 Be but fails to 
describe 1 2 3 4 5 6 He. The zero range force necessary to get agreement in C-Be-Li nuclei has very different parameters 
compared to the parameters fitted by Garrido et al. to reproduce the gap calculated in nuclear matter with the 
finite range Gogny effective interaction. The same Tamm-Dancoff calculation with these two forces shows that 
they are still equivalent in our finite nuclei and gives a good binding energy in 6 He but too weak binding in 
N=8 nuclei. This result is well understood in terms of two-body correlations in the cores. Indeed we know that 
the alpha-particle is well described in Hartree-Fock model. Then two-body correlations in 4 He are unefficient 
while we know for long from shell model calculations that the cores of 12 C and 12 Be and very likely °Li are not 
pure closed shell nuclei as assumed in a pairing model. It is also obvious for 12 Be in our calculation: when it 
is described as a core of 10 Be plus two neutrons its wave function is a mixture of (2s) 2 , (lpi/ 2) 2 and (ld 5 / 2 ) 2 
two-neutron states while in the study of 14 Be it is considered as a pure (2s) 2 or (lpi/ 2) 2 state what yields 
inconsistency of the model. 

The particle-particle RPA model is well adapted to take into account such correlations and indeed the model 
applied to 14 9 C- 11 Li- 12 Be- 14 Be gives now with those realistic forces very good agreement with experimental 
binding energies. It gives also the two-neutron separation energy in the cores. For 12 C- 10 Be and 12 Be (the 
latter being considered as a core in the calculation of 14 Be) the agreement with measurements is also very 
good. However it is too small in 7 Li very likely due to the single proton in the p 3 / 2 shell. The model gives also 
two-neutron and two-neutron hole amplitudes in the wave functions (spectroscopic factors) which are related to 
the amount of two-body correlations introduced in the cores. This effect is always large but larger in 10 Be- 12 Be 
than in 12 C, in qualitative agreement with shell model calculations. 

To get a good two-neutron separation energy in 14 Be and to get a consistent description of 12 Be when it is 
considered as the core of 14 Be or described as 10 Be + two neutrons, one has to assume an inversion of 2s-lp!/ 2 
shells in 12 Be- 13 Be as it is in 10 Li and u Be. This inversion is also suggested by a recent measurement 1331 of 
the transition 2 + (2.1 MeV) —> 0 + (gs) in 12 Be. The B(E2) is found to be the same as for the transition 2 + (3.3 
MeV) —► 0 + (gs) in 10 Bc , suggesting the same effect of particle-phonon couplings in the two systems, n Be- 13 Be, 
therefore the same shells inversion. Furthermore strong particle-particle RPA correlations are known to modify 
the one-neutron mass operator |3^ , |39|| therefore to give further corrections to the one-neutron single energies. 
They may enhance the inversion process studied in refs. [11,34] for 41 Be, even though when adding the two 
contributions due to couplings with phonons and pair vibrations one has to substract the second order term in 
order to eliminate double counting, so that only very mixed RPA states will contribute. From the calculated 
amplitudes one may expect this contribution to be non negligible in 14 Be and 13 Be and even larger in 13 Be than 
in n Be. 

One of us (NVM) is very grateful to P. Schuck for very fruitful discussions which have initiated an important 
part of this work. We also are indebted to G.F. Bertsch for many comments and suggestions during the 
preparation of the manuscript. The second version of the manuscript has been carefully read by F. Bortignon, 
D. Kadi-Hanifi and P. Schuck. We thank them for accepting this tedious job and for their helpful remarks. 


[1] N. Vinli Mau and J.C. Pacheco, Nucl. Phys. A607, 163 (1996). 

[2] G. Audi and A.H. Wapstra, Nucl. Phys. A565, 66 (1993). 

[3] H. Esbensen, G.F. Bertsch and K. Hencken, Phys. Rev. C 56, 3054 (1997) 

[4] X. Campi, J. Martorell and D W.L. Sprung, Phys. Lett. B 41, 443 (1972). 

[5] S. Cohen and D. Kurath, Nucl. Phys. 73, 1 (1965). 

[6] E. Garrido, P. Sarriguren, E. Moya de Guerra and P. Shuck, Phys. Rev. C 60, 064312 (1999). 

[7] J. Decharge and D. Gogny, Phys. Rev. C 21, 1568 (1980). 

[8] J.F. Berger, M. Girod and D. Gogny, Comp. Phys. Comm. 63, 365 (1991). 

[9] G.F. Bertsch and H. Esbensen, Ann. of Phys. (N.Y.) 209 ,327 (1991). 


9 




[10] M. Labiche, F.M. Marques, O. Sorlin and N. Vinh Mau, Phys. Rev. C 60, 027303 (1999). 

[11] N. Vinh Mau, Nucl. Phys. A592, 33 (1995) 

[12] J. Terasaki, H. Flocard, P.H. Heenen and P. Bonche, Nucl. Phys. A621, 706 (1997). 

[13] R. Satchler, L.W. Owen, A.J. Elwyn, G.L. Morgan and R.L. Walter, Nucl. Phys. A112, 1 (1968). 

[14] M. Girod and B. Grammaticos, Phys. Rev. C 27, 2317 (1983). 

[15] T. Suzuki and T. Otsuka, Phys. Rev. C 56, 847 (1997) 

[16] B.A. Brown, International School of Heavy Ion Physics: Exotic Nuclei, edited by R.A. Broglia and P.G. Hansen 
(World Scientific, Singapore, 1998) p.l 

[17] H. Sagawa, T. Suzuki, H. Iwasaki and M. Ishiliara, Phys. Rev. C 63, 034310(2001). 

[18] A. Navin et al., Phys. Rev. Lett. 85, 266 (2000). 

[19] I.J. Thompson and M.V. Zukhov, Phys. Rev. C 53, 708 (1996). 

[20] N. Fukuda, F. Iwamoto and K. Sawada, Phys. Rev. A 135, 932 (1964). 

[21] G. Ripka and R. Padjen, Nucl. Phys. A132, 489 (1969). 

[22] A. Bouyssy and N. Vinh Mau, Phys. Lett. B 35, 269 (1971); Nucl. Phys. A224, 331 (1974). 

[23] S. Pita, These, Universite Paris 6 (2000) 

[24] I.J. Thompson and M.V. Zukhov, Phys. Rev. C 49, 1904 (1994). 

[25] G.F. Bertsch, K. Hencken and H. Esbensen, Phys. Rev. C 57, 1366 (1998) 

[26] R.M. Young et al., Phys. Rev. C 49, 279 (1994). 

[27] T. Kobayashi, Nucl. Phys. A538, 343c (1992). 

[28] M. Zinzer et al., Phys. Rev. Lett. 75, 1719 (1995). 

[29] A.N. Ostrowski et al., Zeit. Phys. Lett. A343, 489 (1992). 

[30] A.V. Belozyorov et al., Nucl. Phys. A636, 419 (1998). 

[31] M. Thoennessen, S. Yokoyama and P.G. Hansen, Phys. Rev. C 63, 014308 (2000). 

[32] F. Ajzenberg-Selove, Nucl. Phys. A506, 1 (1990). 

[33] H. Iwasaki et al., Phys. Lett. B 481, 77 (2000). 

[34] F.M. Nunes et al., Nucl. Phys. A609, 43 (1996). 

[35] P. Descouvemont, Phys. Lett. B 331, 271 (1994). 

[36] P. Descouvemont, Phys. Rev. C 52, 704 (1995). 

[37] Zhongzhou Ren, Gongou Xu, Baoqiu Chen, Zhongyu Ma and W. Mittig Phys. Lett. B 351, 11 (1995). 

[38] N. Vinh Mau, Theory of nuclear structure: Trieste Lectures 1969,(IAEA Vienna) (1970) p.931 

[39] C. Barbieri and W.H. Dickhoff, Phys. Rev. C 63, 034313 (2001). 


10 



